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Aanleiding  van  het  onderzoek 

De  wens  de  specifieke  voordelen  van  glasvezel  versterkte  kunststoffen  op  uitgebreidere  schaal  te 
benutten. 

Doel  van  het  onderzoek 

Het  ontwikkelen  van  kosteneffectieve  oplossingen  voor  ontwerp  en  fabricage  van  composiet 
constracties  in  de  scheepsbouw,  waarmee  deze  constmcties  op  veel  grotere  schaal  toegepast  kunnen 
worden  en  bovendien  een  langere  levensduur  hebben. 


Korte  omschrijving 

Een  methode  uit  de  literatuur  "de  gereduceerde  stijfheidsmethode"  is  bestudeerd  en  gebruikt  om 
van  verschillende  grafiet-epoxy  laminaten  het  begin  van  delaminatie  te  bepalen. 

Conclusies  en  aanbevelingen 

De  gereduceerde  stijfheidsmethode  is  gebruikt  om  het  begin  van  delaminatie  van  een  aantal 
laminaten  van  AS4/3501.6,  T300/934  en  T300/5208  grafiet-  epoxy  platen  te  bepalen.  De  resultaten 
zijn  veelbelovend.  Het  voordeel  is  dat  het  niet  nodig  is  om  gedetailleerde  spanningen  in  de  derde 
richting  en  lokale  spanningen  aan  het  oppervlak  van  de  plaat  te  berekenen.  Ten  einde  de  methode 
bruikbaar  te  maken  voor  de  eindige-elementen  methode  is  het  nodig  situaties  met  varierende 
spanningen  te  onderzoeken.  Ook  is  het  mogelijk  dat  het  begin  van  delaminatie  verschillend  zal  zijn 
voor  de  twee  loodrechte  richtingen  van  het  laminaat.  In  gevallen  waar  beide  (loodrechte) 
spanningen  dicht  bij  de  kritieke  delaminatie  spanning  liggen,  zullen  ze  waarschijnlijk  interactie 
vertonen.  Ook  de  interactie  tussen  matrix  scheuren  en  delaminatie  moet  verder  onderzocht  worden. 
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1  INTRODUCTION 

In  this  report  some  aspects  of  delamination  (separation  of  laminate  plies)  are  considered.  If  a 
laminate  is  loaded  in  its  plane  it  may  happen  that  two  plies  separate.  Because  the  various  plies  in 
the  laminate  contract  differently  perpendicular  to  the  loading,  stresses  in  the  other  directions  arise. 
Near  the  boundary  shear  stresses  are  necessary  to  build  up  these  inplane  stresses.  Moreover  in  the 
third  direction  stresses  arise  to  prevent  the  laminate  to  bend  due  to  the  inplane  stresses.  Due  to  the 
combined  influence  of  shear  forces  and  the  stresses  in  the  third  direction  delamination  may  arise. 
This  may  lead  to  a  significant  loss  of  stiffness.  The  reduced  stiffness  method  has  been  used  to 
determine  the  onset  of  delamination  in  a  laminate.  This  theory  has  the  advantage  that  detailed 
calculations  of  the  stresses  in  the  other  directions  are  not  necessary.  It  will  be  outlined  in  Chapter 
2.  With  the  theory  we  calculate  the  onset  of  delamination  in  a  number  of  different  layups  of  some 
graphite  epoxy  layups,  for  which  data  were  available.  The  results  are  presented  in  Chapter  3. 
Finally  we  discuss  the  results  and  give  some  recommendations  in  Chapter  4. 


2  THEORY 

If  the  plies  of  a  laminate  separate  (called  delamination),  its  longitudinal  laminate  stiffness  (in  the 
loading  direction)  reduces.  The  energy  which  is  released  by  this  stiffness  reduction  is  used  to 
create  the  new  (delamination)  surface.  O’Brien  [1]  used  this  reduced  stiffness  method  to 
determine  the  onset  of  delamination.  The  theory  will  be  explained  using  the  next  figure. 
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/ 

Figure  1  Partially  delaminated  specimen 

A  specimen,  consisting  of  a  number  of  plies  with  different  fibre  directions,  has  been  delaminated  at 
both  sides  over  a  certain  distance.  The  specimen  is  loaded  in  the  longitudinal  (x-)direction 
(perpendicular  to  the  plane  of  this  figure).  The  effective  stiffness  in  this  direction  of  this  specimen 
can  be  approximated  as: 
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EpHEd-Ei)^Ei={Ed-Ei)Dd+Ei 


(1) 


where: 

Ep  :  Stiffness  of  partially  delaminated  specimen  [Pa]; 
£,  :  Stiffness  undamaged  specimen  [Pa]; 

Ed  :  Stiffness  totally  delaminated  specimen  [Pa]; 

2  b  :  Width  of  the  specimen  [m]; 
a  :  Delamination  length  at  each  side  [m]. 

Dd  :  Damage  parameter  for  delamination  [-]. 


This  approximation  is  valid  as  long  as  the  undelaminated  part  is  large  compared  to  the  thickness. 
The  ratio  of  a  and  b  is  a  delamination  damage  parameter  Dd,  which  is  1  for  total  delamination.  It 
will  be  used  for  damage  growth  later  on.  The  amount  of  stored  energy  n  in  the  specimen  is  (the 
strains  in  the  y  and  z-directions  are  unimportant  because  there  is  no  external  load  in  that 
directions): 


n=  fie  C..  e.  dV  =  fieX  dV 
J  2  1  y  j  J  2 

V  V 


(2) 


V  =  2  bit 


where: 

C  :  Stiffness  tensor  [Pa]; 

Ex  :  Strain  in  the  longitudinal  direction  [-]; 

V  :  Volume  [m3]; 

1  :  Length  of  the  specimen  [m]; 

t  :  Thickness  of  the  laminate  [m]. 

The  energy  release  rate  Gd  is: 

r  =_  dn  _  (el Ep bit)  =_  e2xbt  dEp  =  _&l(Ed-E)t 
d  Ha  8(2 at)  2 

A  =2  al 


where: 

A  :  Delaminated  surface  [m2]; 

With  this  equation  the  strain  ec  where  delamination  initiates  can  be  determined  from  the  critical 
energy  release  rate  Gdc  (or  Gdc  from  an  experimental  £c ): 


e.  = 


N 


2  G 


dc 


m  -Ed) 


(4) 


where: 

£c  :  Delamination  initiation  strain  [-]; 

Gdc  :  Energy  release  rate  at  delamination  initiation  [N/m]. 


TNO-Report 

93-CMC-R0786 


Date 

30  August  1996 


Page 

6 


Damage  growth  is  treated  in  the  same  way  as  for  matrix  cracks  as  reported  by  Paas  and  van  den 
Eikhoff  [2].  A  thermodynamic  force  for  delamination  Xd  is  defined  by: 


y  __ 

3dT~ 


3(^X) 


2  A  P'  _  1/77  T?  \  n?  —  d 


3D, 


3D, 


(5) 


where  \|r  is  the  specific  stored  energy. 


If  we  choose  a  dissipation  potential  <j>  linear  in  the  thermodynamic  force: 

v°  - 

A.d  ~ - 

t 


(6) 


with: 

Ad  :  Material/constmction  parameter  for  delamination  expansion  [-]; 
Xd°  :  Threshold  value  for  delamination  growth  [N/m2]; 


This  potential  can  be  used  to  determine  the  delamination  damage  as  a  function  of  strain  in  an 
analog  way  as  in  [2]: 


GJ 


Dd  =  ^-=AA-1-1)=AA -JL  -1)-A,(4-1) 


''  b  ‘“'X! 


(7) 


'dc 


The  "strain  at  total  delamination  £d"  can  be  determined  by  substituting  Dd=  1  (a=b)  in  (7): 


ed=e- 


N 


l+± 


(S) 


The  values  for  E{  en  Ed  have  been  calculated  for  a  number  of  layups  with  PLAMOR.  During  this 
project  it  became  clear  that  there  is  a  lot  of  confusion  which  stiffness  term  has  to  be  used.  We  will 
now  shortly  indicate  this  problem. 


For  a  general  laminate  we  have  the  following  relation  between  laminate  forces  and  moments  on  the 
one  hand  and  the  midplane  strain  and  curvature  at  the  other  hand: 


N 

A 

B 

£ 

M 

BT 

D 

K 

where: 

N :  Forces  on  the  laminate  [A^,  Ny,  AyT; 

M:  Moments  on  the  laminate;  [Afx,  My,  M^] T; 

e:  Strain  in  the  midplane;  [ex,  £y,  £xy]T; 

k:  Curvature  of  the  midplane  [k^  Ky,  k^,]7;. 
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For  the  reduced  stiffness  we  are  only  interested  in  the  reaction  force  in  the  load  direction  on  a 
strain  in  that  direction.  However  this  reaction  force  will  depend  on  the  other  boundary  conditions 
of  the  laminate.  If  we  constrain  ourselves  to  symmetrical  laminates  which  are  partly  delaminated 
from  the  edges  and  loaded  in  one  direction,  than  we  have  the  following  boundary  conditions: 

A)  Natural  boundary  conditions  in  the  delaminated  part 

The  delaminated  parts  can  freely  extend  in  the  y-direction  (in  the  plane  of  the  laminate 
perpendicular  to  the  load)  (Ny=0)  and  they  can  bend  freely  about  the  x-axis  {My- 0).  It  can 
not  shear  around  the  z-axis  (exy=0)  and  probably  there  will  be  no  twist  (1^=0),  furthermore 
the  faces  will  remain  approximately  straight  (^=0).  The  details  of  the  determination  of  the 
longitudinal  (x)  laminate  stiffness  are  given  in  Appendix  1. 

B)  Alternating  delamination  boundary  conditions 

From  experiments  it  follows  that  sometimes  the  delamination  switches  up  and  down 
between  two  symmetric  interfaces  in  the  length  of  the  specimen.  O’Brien  [1]  believes  that 
"Delaminations  grow  in  this  matter  to  reduce  the  effect  of  bending-extension  coupling". 
This  means  that  internally  there  may  be  a  (opposite)  curvature  in  two  adjacent  parts  of  the 
laminate,  giving  a  nett  effect  of  approximately  no  curvature.  The  longitudinal  laminate 
stiffness  is  given  in  Appendix  2. 

C)  Free  boundary  conditions 

Schellekens  [3]  uses  a  formula  for  the  effective  stiffness  which  is  equivalent  with  the 
assumption  that  all  the  moments  and  Ny  and  are  zero.  From  (9)  the  curvatures  can  than 
be  eliminated  giving  the  equations  of  Appendix  3. 

D)  Curvature  free  boundary  conditions 

O’Brien  [1]  found  that  using  the  free  boundary  conditions  gives  results  that  are  not  stiff 
enough  and  he  proposes  to  discard  the  bending  extension  coupling,  which  is  the  same  as 
demanding  that  all  the  curvatures  are  zero,  see  Appendix  4. 

It  is  the  opinion  of  the  author  that  the  boundary  conditions  A)  and  B)  are  the  right  ones  to  use, 
whichever  gives  the  lowest  delamination  initiation  strain. 


3  EXAMPLES 

In  this  chapter  we  apply  the  theory  of  the  previous  chapter  for  a  number  of  different  graphite 
epoxy  laminates  to  calculate  the  delamination  initiation  strain  and  for  which  experimental  or 
numerical  data  were  available.  The  program  PLAMOR  developed  at  the  TU  Delft  has  been  used  to 
determine  the  stiffnesses  of  the  laminates.  In  Section  3.1  we  focus  on  the  comparison  with 
theoretical  calculations.  Delamination  growth  is  the  major  topic  in  Section  3.2.  Thickness  effects  is 
the  topic  in  Section  3.3.  Sections  3.4  and  3.5  show  that  interaction  with  matrix  cracking  or  fibre 
break  is  important  for  certain  laminates. 


3.1  Comparison  with  theory  for  [±25/90]s  layup  of  AS4/3501.6  graphite  epoxy 

In  this  section  we  compare  the  reduced  stiffness  results  with  results  of  calculations  by  Schellekens 
en  de  Borst  [3]  with  a  [±25/90]s  layup  of  AS4/3501.6  graphite  epoxy.  They  gave  the  following 
material  data: 

Ei 

e2 


:  140  GPa; 
:  11  GPa; 
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G12  :  5.5  GPa; 

vn  :  0.29; 

ply  thickness  :  0.132  mm. 

For  a  number  of  layups  the  effective  Ex  has  been  determined  with  PLAMOR  and  the  formulas  of 

Appendix  1  and  3.  The  results  are  given  in  the  next  table  for  Free  boundary  conditions  (App.  3) 
used  in  [3]  and  "Natural  boundary  conditions"  (App  1),  which  should  have  been  used  according  to 
the  author. 


It  is  obvious  that  the  boundary  conditions  have  an  enormous  influence  on  the  longitudinal  laminate 
stiffness,  the  free  conditions  being  much  less  stiff.  By  taking  combinations  of  two  layups  (e.g. 
[25/-25/90/90]  and  [25/-25])  and  the  corresponding  value  for  E„  Ed  for  a  complete  delaminated 
layup  can  be  determined  with  delaminations  at  different  interfaces: 


GO) 


where: 

Exn  :  Ex  of  the  part  with  n  plies; 

E”  :  Ex  of  the  part  with  m  plies; 


Next  we  can  calculate  with  (4)  the  critical  strain  where  delamination  would  start  in  the  different 
interfaces.  The  results  are  given  in  the  next  table.  We  used  Gdc= 350  N/m  [3]  (supposed  to  be 
independent  of  the  angle  between  the  fibres  in  adjacent  plies). 


Boundary  condition 

Free*5 

Natural  *) 

Layup 

Ed  [GPa] 

ec  H 

Ed  [GPa] 

M-l 

[25/-25/90/90/-25]+[25] 

33.6 

.0052 

61.6 

.00141 

[25/-25/90/90M-25/25] 

30.8 

.0050 

55.1 

.00900 

[25/-25/90] + [90/-25/25 ] 

36.1 

.0054 

55.5 

.00917 

[25/-25/90/90]+[-25/25] 

Alternating  delam*); 

54.9 

.00891 

See  Chapter  2  and  the  Appendices 
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From  this  table  follows  that  there  is  little  preference  for  delamination  at  the  different  interfaces  (£c 
is  almost  equal),  but  delamination  between  a  90-layer  and  a  -25-layer  will  probably  start  first, 
preferably  alternating  between  the  two  symmetric  interfaces.  The  influence  of  the  boundary 
conditions  is  tremendous,  approximately  a  factor  2  in  the  delamination  initiation  strain. 

In  [3]  it  was  found  that  for  a  (non-linear)  fracture  mechanics  approach  with  plane  strain  elements 
and  3-D  line  interface  elements  a  critical  strain  was  found  of  approximately  .009  for  a 
delamination  in  the  midplane.  For  the  natural  boundary  conditions  and  the  reduced  stiffness 
method  we  find  almost  the  same  initiation  strain.  In  order  to  get  the  same  answer  with  (4)  and  the 
free  boundary  conditions  one  should  use  a  Gdc  of  970  N/m.  Schellekens,  who  used  these 
conditions  for  the  reduced  stiffness,  concluded  that  "this  cast  some  doubt  on  the  merit  of  the 
reduced  stiffness  approach  in  predicting  the  onset  of  delamination".  We  now  can  conclude  that  the 
results  of  this  method  and  the  Finite  Element  Method  match  perfectly,  provided  the  same  boundary 
conditions  are  used  in  both.  We  must  also  conclude  that  it  is  dangerous  to  rely  on  symmetry 
conditions  because  unlike  the  assumptions  in  the  Finite  Element  calculations  the  delamination  will 
probably  not  start  in  the  midplane  of  the  laminate  and  may  not  even  be  on  the  same  side  of  the 
midplane  over  the  whole  length  of  the  specimen. 


3.2  Delamination  growth  in  [±30±30/90,,/JS  layups  of  T300/5208  graphite  epoxy 

Now  we  deal  with  an  example  of  [±30±30/90,,Js  layups  of  T300/5208  graphite  epoxy  as  presented 
by  O’Brien  [1].  The  material  data  were  given  as: 

Et  :  138  GPa; 

E2  :  15  GPa; 

G12  :  5.9  GPa; 

v12  :  0.21; 

ply  thickness  :  0.137  mm. 

From  this  data  we  calculated,  for  the  natural  boundary  conditions: 


Layup 

Ex  [GPa] 

[±30]2 

53.8 

[903] 

15.0 

[±30/±30/903] 

41.6 

[±30/±30/90m] 

47.0 

[±30/±30/90s/±30/±30] 

58.2 
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and: 


Layup 

E, 

[GPa] 

Ed 

[GPa] 

£c  H 

ec  [-] 

experim. 

[±30/±30/903]+[±30]2 

58.2 

46 

.00376 

.0035 

[±30]2+[903]+[±30]2 

58.2 

43.2 

.00479  * 

.0035 

[±30/±30/901,/J+[901,/y±30/±30] 

58.2 

47.0 

.00392 

.0035 

[±30/±30/903]+[±30]2  (alternating) 

58.2 

45.9 

.00374 

.0035 

Here  two  delamination  surfaces  arise. 


A  value  of  Gdc=  130  N/m  has  been  determined  in  [1]  to  fit  the  experimental  delamination  initiation 
and  has  been  used  with  (4)  to  determine  the  delamination  initiation  strain.  This  value  is  much 
smaller  than  the  value  used  in  the  previous  examples  for  AS4/3501.6. 

A  value  for  G&,  independent  of  the  experimental  data,  can  be  determined  as  follows: 

The  critical  energy  release  rate  Glc  for  transverse  matrix  cracking  was  found  to  be  [2]: 

□  64  N/m  for  AS4-3501.6 

□  31  N/m  for  T300/934  . 

Assuming  that  the  same  ratio  will  hold  for  the  critical  release  rate  in  delamination,  Gdc  =  170  N/m 
(31/64*350)  can  been  used,  giving  an  (alternating)  delamination  initiation  strain  of  .0043  in 
reasonable  accordance  with  the  experimental  data. 


[1]  also  gives  the  development  of  Gd  versus  the  relative  delaminated  surface  (i.e.  the  delamination 
damage)  as  shown  in  the  next  Figure. 
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T  Gd  [N/m] 


relative  delaminated  surface  Dd  [-] 

Figure  2  Delamination  resistance  curve  for  [±30/±30/903/ ±30/ ±30]  T300-5208  graphite- 
epoxy  laminates. 


Using  Equation  (7)  we  see  that  a  large  part  of  the  resistance  curve  (up  to  a  damage  of  0.7)  can  be 
fitted  using  Ad  =  0.5  and  Gdc  =  160  N/m  (  which  is  very  close  to  the  value  of  170  N/m  which  was 
found  as  a  prediction).  Gd  increases  during  delamination  growth.  From  this  figure  it  can  also  be 
seen  that  delamination  growth  starts  in  the  experiments  for  Gd  between  130  and  200  N/m.  For 
damages  above  0.7  the  resistance  curves  bend  upwards. 


3.3  Thickness  influence  in  [±45J0J90JS  layups  of  T300/5208  graphite  epoxy 

O’Brien  [1]  also  gives  results  for  [±45n/0n/90n]s  layups  of  T300/5208  graphite  epoxy.  The  material 
data  were  the  same  as  in  the  previous  example,  however  the  ply  thickness  was  0.15  mm. 
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For  the  various  layups  we  find: 


Layup 

Ex  [GPa] 

[45] 

20.6 

[±45] 

20.9 

[±45/0] 

60.4 

[±45/0/90] 

49.3 

[±45/0/902] 

43.0 

[±45/0/902/0] 

59.4 

[±45/0/902/0/-45] 

58.7 

[  ±4 5/0/ 902/ 0/±45  ] 

55.7 

We  determined  the  critical  strain  in  the  various  interfaces  for  the  same  Gdc=130  N/m: 


Layup 

E;  [GPa] 

Ed  [GPa] 

M-] 

[45]+[±45/0/902/0/-45] 

55.7 

53.9 

.0108 

[±45]+[±45/0/902/0] 

55.7 

49.8 

.0060 

[±45/0]+[±45/0/902] 

55.7 

49.5 

.0059 

[±45/0/90] +[±45/0/90] 

55.7 

49.3 

.0058 

From  this  table  it  follows  that  delamination  most  likely  will  start  in  the  surface  of  symmetry.  For 
this  layup  an  alternating  delamination  was  not  predicted.  For  three  different  thicknesses  we 
compare  the  results  with  the  experimental  delamination  initiation  strains  given  in  [1]. 


n 

ec  [-] 

M-] 

theory 

experiment 

i 

.0058 

.0058 

2 

.0041 

.0047 

3 

.0033 

.0036 

We  can  conclude  that  the  correspondence  between  experimental  value  and  predicted  value  is  good 
although  the  Gic  from  another  layup  has  been  used. 
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3.4  Damage  interaction  in  [±25/90 Js  layups  of  T300/934  graphite  epoxy 

Crossman  and  Wang  [4]  give  results  for  [±25/90n]s  layups  of  T300/934  graphite  epoxy  for 
different  thicknesses  of  the  90-layer.  We  used  their  material  data: 

£,  :  145  GPa; 

E2  :  11.7  GPa; 

G12  :  6.5  GPa; 

v12  : 0.3. 

No  thickness  was  reported,  we  used  t=0.132  mm. 

For  a  number  of  layups  the  effective  Ex  is  calculated  with  PLAMOR  for  the  natural  boundary 
conditions  and  for  the  alternating  boundary  conditions.  From  these  values  the  delamination 
initiation  strain  can  be  calculated  and,  with  the  initial  stiffness,  the  delamination  initiation  stress. 
The  results  for  an  alternating  delamination  are  presented  in  the  next  table,  because  the  initial 
stresses  are  the  lowest.  The  results  are  compared  with  the  experimental  stresses  in  the  next  table 
using  the  same  G^BO  N/m  as  in  the  previous  example  (but  a  different  epoxy).  The  layup  with 
the  lowest  initiation  stress  has  been  presented.  It  appears  that  delamination  is  predicted  at  the  90/25 
interface  for  all  n.  For  n< 3  delamination  in  the  surface  of  symmetry  is  almost  as  likely.  This  is  in 
agreement  with  the  experiments  where  both  arise  (sometimes  even  alternating)  for  n< 3.  For  n> 3 
two  alternating  delaminations  arise  are  seen  in  the  experiments.  For  those  thick  90-layers  delami¬ 
nation  starts  along  the  transverse  cracks  indicating  that  these  two  damage  modes  interact. 

From  the  calculations  follows,  that  there  is  little  difference  between  delamination  at  the  symmetry 
plane  and  at  the  -25/90  plane. 


Layup 

Et  [GPa] 

Ed  [GPa] 

cc  [MPa] 

oc  [MPa] 
experim. 

[±25/90]+[±25] 

77.6 

65.9 

449 

425-457 

[±25/902]+[±25] 

69.0 

58.3 

383 

409-413 

[±25/90J+[±25] 

55.8 

48.1 

316 

311-332 

[±25/906]+[±25] 

47.3 

41.3 

271 

258-287 

[±25/908]+[±25] 

41.5 

36.6 

239 

200-216 

[±25/9012]+[±25] 

34.2 

30.5 

197 

125-135 

[±25/9016]+[±25] 

29.8 

26.8 

172 

92.3-97.8 

[±252/904)+[±252] 

69.0 

58.3 

272 

296-308 
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The  calculated  initial  stiffness  is  compared  with  the  experimental  stiffness  in  the  next  Figure  3. 
The  correspondence  is  very  good. 


Delamination  T300/934 


[+-25/90n]s 


□  Exper . 


number  of  90-plies  n 
+  Calc. 


Figure  3  Tangent  modulus  for  [±25/90n]s  T300/934. 


The  dependence  of  the  delamination  initiation  stress  with  the  number  of  90-plies  is  shown  in 
Figure  4  and  is  compared  with  the  experimental  range  (lower-upper).  The  general  trend  is  the  same 
for  experiment  and  calculation,  but  the  mean  of  the  calculated  stresses  is  16%  higher  than  the 
mean  of  the  experimental  values.  Especially  for  thick  laminates  the  experimental  values  are  lower, 
probably  because  of  the  interaction  with  the  transverse  matrix  cracks. 

Delamination  initiation  strain  versus  the  number  of  the  90-plies  is  shown  in  Figure  5  together  with 
the  experimental  range.  The  predicted  values  are  very  good  up  to  n=4.  It  can  be  seen  that  the 
experimental  initiation  strain  drops  for  thick  layers,  whereas  the  calculated  value  remains  almost 
constant.  For  small  n  (up  to  4)  the  (experimental)  initiation  strain  is  7%  below  the  ultimate  strain 
of  a  [±25]s  laminate;  for  n=5  and  6  it  is  7%  below  the  ultimate  strain  of  a  [90] 16  laminate. 


eps  [mm/m]  slgmo  CMPa] 
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Figure  5 


Del  ami  nation  T300/934 


Delamination  stress  of  [±25/90n]s  T300/934. 
Lower-Upper:  Experimental  range. 


Del  ami  nation  T300/934 


o  2  4  E 

number  of  90-piles  n 

□  Lower  +  Upper  O  Calc. 

Delamination  strain  of  [±25/90n]s  T300/934. 
Lower-Upper:  Experimental  range. 
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This  drop  in  initiation  strain  is  probably  due  to  the  interaction  with  the  transverse  cracks.  Because 
of  these  cracks  the  stiffness  of  the  laminate  drops  considerably,  especially  for  thick  90-layers, 
giving  larger  strains  than  the  ones  following  from  the  original  stiffness.  This  can  also  explain  why 
the  results  for  the  delamination  initiation  stress  agree  better,  because  the  initial  stiffness  has  been 
used  to  calculate  the  initial  strains  from  the  experiments.  Also  the  strain  in  the  outer  layers  is 
concentrated  somewhat  over  the  transverse  crack  giving  a  faster  initiation. 

With  the  method  outlined  in  [2]  we  also  calculated  the  onset  of  transverse  matrix  cracking.  The 
results  are  presented  in  the  next  table. 


Layup 

oc  [MPa] 
transv.  crack 

ac  [MPa] 
Experimental 

[25/-25/90Js 

250 

207-243 

[25/-25/903]s 

171 

147-167 

[25/-25/904]s 

131 

127 

[25/-25/906]s 

88 

98-120 

[25/-25/908]s 

66 

85.4-97.8 

[252/-252/902]s 

320 

296-308 

The  correspondence  is  good,  although  the  differences  in  critical  stresses  are  somewhat  smaller  in 
the  experimental  data  then  in  the  calculated  values.  For  the  thin  specimens  the  onset  of  transverse 
matrix  cracking  is  for  slightly  lower  strains,  so  when  delaminations  starts  there  are  only  a  few 
cracks.  For  thick  laminates  however  transverse  matrix  cracking  starts  at  much  lower  strains.  This 
means  that  there  are  much  more  cracks,  which  will  influence  the  stiffness  and  the  delamination 
initiation. 


3.5  Damage  interaction  in  [0/90  Js  layups  of  T300/934  graphite  epoxy 

Wang,  Kishore  and  Li  [5]  give  experimental  results  for  [02/90Js  crossply  layups.  The  material 
data  were  as  given  in  the  previous  example.  An  estimation  for  the  delamination  initiation  strain  has 
been  calculated.  The  results  are  compared  with  the  experimental  results  in  the  next  table.  Again  we 
used  Gac  of  130  N/m.  The  ply  thickness  was  reported  to  be  0.132  mm. 


Layup 

Ej  [GPa] 

Ed  [GPa] 

M-] 

Ec  [mm/m] 
experim. 

[0/0/90]+[90/0/0] 

101 

100.7 

.030 

>.0123 

[0/0/903]+[90/0/0]  11 

78.7 

78.5 

.035 

>.0100 

[0/0/907]+[90/0/0]  *> 

56.4 

56.3 

.030 

.0080 

[0/0] +[0/0/90 , 6] + [0/0] 2) 

38.5 

38.5 

.038 

.0041 

1}  two  alternating  delaminations  2)  two  continuous  delaminations 
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From  the  calculations  follows  that  delamination  is  most  likely  on  the  plane  of  symmetry  for 
[02/90]s  and  near  the  0/90  interface  for  the  thicker  specimens.  For  the  thickest  laminate  two 
continuous  delaminations  at  the  0/90  interface  are  predicted.  The  order  of  the  predicted  initiation 
strain  is  much  higher  than  the  strains  where  matrix  cracking,  0-ply  splitting  and  fibre  break  occur. 
The  calculated  values  are  nearly  the  same  for  all  thicknesses.  It  is  likely  that  in  this  case  there  is  a 
strong  interaction  with  the  other  damage  modes. 


4  DISCUSSION  AND  RECOMMENDATIONS 

In  the  previous  chapter  we  have  applied  the  reduced  stiffness  method  on  layups  of  three  different 
graphite  epoxies.  The  results  are  very  good  and  the  advantage  is  that  it  is  not  necessary  to  calculate 
detailed  stresses  in  the  other  directions  and  to  calculate  local  stresses  at  the  side  surfaces  of  the 
plate. 

We  found  that  the  discrepancy  between  Finite  Element  calculations  and  the  reduced  stiffness 
method  observed  by  Schellekens  was  caused  by  the  fact  that  he  used  different  boundary  conditions 
in  both  approaches.  If  the  correct  (natural,  see  Appendix  1)  boundary  conditions  are  used  the 
results  match  very  good. 

From  the  calculations  also  follows  that  with  the  same  energy  release  rate  Gdc  the  results  for 
different  lay-ups  can  be  predicted.  Only  for  situations  where  the  delamination  initiation  strain  is 
higher  than  the  strain  where  other  damages  like  matrix  cracking  and  fibre  cracking  initiate,  the 
results  are  poor  because  there  is  an  interaction  with  these  other  damages. 

The  reduced  stiffness  method  can  very  well  be  used  in  design.lt  is  not  yet  clear  whether  it  is 
possible  to  apply  the  method  in  a  rigorous  way  in  the  Finite  Element  Method  in  situations  with 
varying  stresses. 

It  is  possible,  that  the  onset  of  delamination  will  be  different  for  the  stresses  in  the  two 
perpendicular  directions  of  the  laminate.  In  cases  with  one  large  stress  component  there  is  probably 
no  problem,  but  in  cases  where  both  stresses  are  near  to  the  delamination  initiation  stress,  the  two 
will  interact.  This  problem  of  interaction  has  to  be  solved  before  an  implementation  in  a  Finite 
Element  Method  can  be  used  in  a  rigorous  way. 
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Appendix  1  Longitudinal  laminate  stiffness  for 

"Natural  boundary  conditions" 


For  the  "natural  boundary  conditions"  discussed  in  Chapter  2  we  can  derive  the  longitudinal 
laminate  stiffness.  First  rearrange  the  relation  between  forces/moments  and  strains/curvatures  such 
that  we  have  second  the  terms  with  zero  forces/moments  and  last  the  term  with  zero 
strains/curvatures : 


0 

0 

N 

K 

MV. 


e  s T  s T 

J12  ^13 


S'  S'  c 7 

u12  u22  J23 


^13  ^23  ^33 


Because  we  are  only  interested  in  Nx  and  ex  we  eliminate  £y  and  K*,  giving: 


T 


where: 

T:  Total  laminate  thickness. 

This  equation  gives  the  required  longitudinal  laminate  stiffness  term. 


Appendix  2  Longitudinal  laminate  stiffness  for 

"Alternating  delamination  boundary  conditions" 

For  the  "alternating  delamination  boundary  conditions"  discussed  in  Chapter  2  we  can  derive  the 
longitudinal  laminate  stiffness.  Again  rearrange  the  relation  between  forces/moments  and 
strains/curvatures  such  that  we  have  second  the  terms  with  zero  forces/moments  and  last  the  term 
with  zero  strains/curvatures: 


r\- 

r 

"] 

~  _“1 

* 

ST  ST 

°12  °13 

M 

K 

L  XJ 

L 

0 

s" 

L° . 

Sn 

S  c T 

*z2  J23 

K 

L  yj 

_ l 

s* 

523  S33 

0 

_0J_ 

Again  we  eliminate  £y  and  K*,  giving: 
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In  the  delaminated  laminate  we  now  have  the  situation  as  indicated  in  the  following  figure  of  a 
deformed  laminate: 


We  have  two  parts  of  the  delaminated  laminate  on  top  of  each  other.  We  assume  that  the  surface 
between  the  two  opposite  delaminated  areas  remains  straight.  In  that  case  we  get  the  following 
relation  between  the  local  midplane  strains  and  the  global  midplane  strain  and  curvature: 

T2 

e  =  e  -  k  _ 

1  2 

Ty 

£,  =  8  +  K  __ 

2  2 

where: 

8:  strain  at  the  midplane  of  the  total  laminate; 

k:  curvature  in  the  delaminated  parts; 

8;:  strain  in  the  midplane  of  part  i; 

T{.  Thickness  of  part  i. 

Now  the  forces  and  moments  on  the  two  parts  can  be  calculated: 

"  Em  e  +  -  y  E'} 

N2  -  e  +  <S*  +  -y  E2) 

e  +  V 


It  is  now  possible  to  have  no  external  moment  and  a  mean  curvature  over  adjacent  parts  of  zero. 
This  means  that: 
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T  T 
M  +  M  =  N.  -1  -  N,  _i 
1  2  1  2  2  2 


We  now  eliminate  the  curvature: 

K 

h 

Sa\ 

Svi 


si 

Sa2  +  2  S41  +  f(s\ 


from  which  the  longitudinal  laminate  stiffness  follows. 
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Appendix  3  Longitudinal  laminate  stiffness  for 

"moment  free  boundary  conditions" 

For  the  "moment  free  boundary  conditions  discussed  in  Chapter  2  we  can  derive  the  longitudinal 
laminate  stiffness. This  is  done  by  eliminating  the  curvatures  from  (9)  giving: 

N  =  [A  -  Bt  D-1  b]z 

from  which  the  strain  can  be  calculated  with  matrix  inversion: 

E  =  [A  -  BT  D B]'1  N  =  A'  N 

With  N_xy= 0  the  longitudinal  laminate  stiffness  follows: 

1 


o 
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Appendix  4  Longitudinal  laminate  stiffness  for 

’’curvature  free  boundary  conditions" 

For  the  "curvature  free  boundary  conditions  discussed  in  Chapter  2  we  can  derive  the  longitudinal 
laminate  stiffness.  This  is  done  by  matrix  inversion  using  (9)  giving: 

£=  A'1  N 


The  longitudinal  laminate  stiffness  follows  with  A^=0: 
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